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Abstract. In the timescape scenario cosmic acceleration is understand as an apparent effect, due to gravitational energy 
gradients that grow when spatial curvature gradients become significant with the nonlinear growth of cosmic structure. 
This affects the calibratation of local geometry to the solutions of the volume-average evolution equations corrected by 
backreaction. In this paper I discuss recent work on defining observational tests for average geometric quantities which can 
distinguish the timescape model from a cosmological constant or other models of dark energy. 

Keywords: dark energy, theoretical cosmology, observational cosmology 
PACS: 98.80.-k 98.80.Es 95.36,+x 98.80Jk 



INTRODUCTION 

I will discuss some recent results on observational tests yl of a model cosmology, which represents a new approach 
to understanding the phenomenology of dark energy as a consequence of the effect of the growth of inhomogeneous 
structures. The basic idea, outlined in a nontechnical manner in ref. [2], is that as inhomogeneities grow one must 
consider not only their backreaction on average cosmic evolution - as discussed by other contributors to this volume 
- but also the variance in the geometry as it affects the calibration of clocks and rods of ideal observers. Dark energy 
is then effectively realised as a misidentification of gravitational energy gradients. 

Although the standard Lambda Cold Dark Matter (ACDM) model provides a good fit to many tests, there are 
tensions between some tests, and also a number of puzzles and anomalies. Furthermore, at the present epoch the 
observed universe is only statistically homogeneous once one samples on scales of 150-300 Mpc. Below such scales 
it displays a web-like structure, dominated in volume by voids. Some 40%-50% of the volume of the present epoch 
universe is in voids with dp/p ~ —1 on scales of 30/i _1 Mpc JH], where h is the dimensionless parameter related to the 
Hubble constant by H Q = 100/ikm sec -1 Mpc -1 . Once one also accounts for numerous minivoids, and perhaps also a 
few larger voids, then it appears that the present epoch universe is void-dominated. Clusters of galaxies are spread in 
sheets that surround these voids, and in thin filaments that thread them. 

One particular consequence of a matter distribution that is only statistically homogeneous, rather than exactly 
homogeneous, is that when the Einstein equations are averaged they do not evolve as a smooth Friedmann-Lemaitre- 
Robertson- Walker (FLRW) geometry. Instead the Friedmann equations are supplemented by additional backreaction 
terms |4J. Whether or not one can fully explain the expansion history of the universe as a consequence of the growth 
of inhomogeneities and backreaction, without a fluid-like dark energy, is the subject of ongoing debate (5J. 

Elsewhere in this volume, Peebles ||6|] provides some of the arguments that have been presented against backreaction. 
His line of reasoning is that of a plausibility argument: if we assume a FLRW geometry with small perturbations, and 
estimate the magnitude of the perturbations from the typical rotational and peculiar velocities of galaxies, then the 
corrections of inhomogeneities are consistently small. This would be a powerful argument, were it not for the fact that 
at the present epoch galaxies are not homogeneously distributed. The Hubble Deep Field reveals that galaxies were 
close to being homogeneous distributed at early epochs, but following the growth voids at redshifts z < 1 that is no 
longer the case today. Therefore galaxies cannot be consistently treated as randomly distributed gas particles on the 
30h~ l Mpc scales y|] that dominate present cosmic structure below the scale of statistical homogeneity. 

Over the past few years I have developed a new physical interpretation of cosmological solutions within the Buchert 
averaging scheme O, & HI] • I start by noting that in the presence of strong spatial curvature gradients, not only 
should the average evolution equations be replaced by equations with terms involving backreaction, but the physical 
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interpretation of average quantities must also account for the differences between the local geometry and the average 
geometry. In other words, geometric variance can be just as important as geometric averaging when it comes to the 
physical interpretation of the expansion history of the universe. 

I proceed from the fact that structure formation provides a natural division of scales in the observed universe. As 
observers in galaxies, we and the objects we observe in other galaxies are necessarily in bound structures, which 
formed from density perturbations that were greater than critical density. If we consider the evidence of the large 
scale structure surveys on the other hand, then the average location by volume in the present epoch universe is in a 
void, which is negatively curved. We can expect systematic differences in spatial curvature between the average mass 
environment, in bound structures, and the volume-average environment, in voids. 

Spatial curvature gradients will in general give rise to gravitational energy gradients. Physically this can be un- 
derstood in terms of a relative deceleration of expanding regions of different densities. Those in the denser region 
decelerate more and age less. Since we are dealing with weak fields the relative deceleration of the background is 
small. Nonetheless even if the relative deceleration is typically of order 10~ 10 ms~ 2 , cumulatively over the age of the 
universe it leads to significant clock rate variances H. I proceed from an ansatz that the variance in gravitational en- 
ergy is correlated with the average spatial curvature in such a way as to implicitly solve the Sandage-de Vaucouleurs 
paradox that a statistically quiet, broadly isotropic, Hubble flow is observed deep below the scale of statistical ho- 
mogeneity. In particular, galaxy peculiar velocities have a small magnitude with respect to a local regional volume 
expansion. Expanding regions of different densities are patched together so that the regionally measured expansion 
remains uniform. Such regional expansion refers to the variation of the regional proper length, t r = Y 1 / 3 , with respect 
to proper time of isotropic observers (those who see an isotropic mean CMB). Although voids open up faster, so that 
their proper volume increases more quickly, on account of gravitational energy gradients the local clocks will also tick 
faster in a compensating manner. 

Details of the fitting of local observables to average quantities for solutions to the Buchert formalism are described 
in detail in refs. QSl- Negatively curved voids, and spatially flat expanding wall regions within which galaxy clusters 
are located, are combined in a Buchert average 

/v(0+./w(r) = l, (1) 

where f w (t) = / w ia w 3 /a 3 is the wall volume fraction and f v (t) — / v i«v 3 /a 3 is the void volume fraction, Y = Ycv' 
being the present horizon volume, and f wl , f vl and f{ initial values at last scattering. The time parameter, t, is the 
volume-average time parameter of the Buchert formalism, but does not coincide with that of local measurements in 
galaxies. In trying to fit a FLRW solution to the universe we attempt to match our local spatially flat wall geometry 

ds 2 = -dT 2 + fl w 2 (T) [dTJ 2 +TJ 2 dn 2 ] . (2) 

to the whole universe, when in reality the rods and clocks of ideal isotropic observers vary with gradients in spatial 
curvature and gravitational energy. By conformally matching radial null geodesies with those of the Buchert average 
solutions, the geometry (O may be extended to cosmological scales as the dressed geometry 

ds 2 = -dT 2 + a 2 (T) [di) 2 + r 2 w (f],T)da 2 ] (3) 

where a = y~ 1 a, y — $z is the relative lapse function between wall clocks and volume-average ones, df] — dt / a = dr / a, 

and r w = y(l -/ v ) 1/3 /wi _1/3 i?w(i?, t), where 77„ is given by integrating d?7 w = f wi 1/3 di) /[f(l - f v ) l/3 } along null 
geodesies. 

In addition to the bare cosmological parameters which describe the Buchert equations, one obtains dressed pa- 
rameters relative to the geometry (f3j). For example, the dressed matter density parameter is Q. M = y 3 ^^, where 

£2 M = SnGp a 3 / (3H 2 a 3 ) is the bare matter density parameter. The dressed parameters take numerical values close 
to the ones inferred in standard FLRW models. 



APPARENT ACCELERATION AND HUBBLE FLOW VARIANCE 

The gradient in gravitational energy and cumulative differences of clock rates between wall observers and volume 
average ones has important physical consequences. Using the exact solution obtained in ref. [80, one finds that a 
volume average observer would infer an effective deceleration parameter q = —a/(H 2 a) =2(1 — f v ) 2 /(2 + / v ) 2 , 



which is always positive since there is no global acceleration. However, a wall observer infers a dressed deceleration 
parameter 

1 d 2 a _ - (1 - fy) (8/ v 3 + 39 f v 2 - 12/ v - 8) 
9 H2a ^ (4 + /v + 4/ v 2 ) 2 ' ^ 

where the dressed Hubble parameter is given by 

H = a- 1 ^a = fH-f = fH-r 1 i- T ?- (5) 

At early times when f v — > the dressed and bare deceleration parameter both take the Einstein-de Sitter value 
q ~q ~ j. However, unlike the bare parameter which monotonically decreases to zero, the dressed parameter becomes 
negative when / v ~ 0.59 and q — > - at late times. For the best-fit parameters [ 10] the apparent acceleration begins at 
a redshift z — 0.9. 

Cosmic acceleration is thus revealed as an apparent effect which arises due to the cumulative clock rate variance 
of wall observers relative to volume-average observers. It becomes significant only when the voids begin to dominate 
the universe by volume. Since the epoch of onset of apparent acceleration is directly related to the void fraction, f v , 
this solves one cosmic coincidence problem. 

In addition to apparent cosmic acceleration, a second important apparent effect will arise if one considers scales 
below that of statistical homogeneity. By any one set of clocks it will appear that voids expand faster than wall regions. 
Thus a wall observer will see galaxies on the far side of a dominant void of diameter 30h~ 1 Mpc recede at a rate greater 
than the dressed global average H Q , while galaxies within an ideal wall will recede at a rate less than H Q . Since the 
uniform bare rate H would also be the local value within an ideal wall, eq. © gives a measure of the variance in the 
apparent Hubble flow. The best-fit parameters IllOll give a dressed Hubble constant H = 61.7+J'jkm sec -1 Mpc -1 , 

and a bare Hubble constant H Q = 48.2^2 4km sec -1 Mpc -1 . The present epoch variance is 17-22%. 

Since voids dominate the universe by volume at the present epoch, any observer in a galaxy in a typical wall region 
will measure locally higher values of the Hubble constant, with peak values of order 72km sec -1 Mpc -1 at the 30/z -1 
Mpc scale of the dominant voids. Over larger distances, as the line of sight intersects more walls as well as voids, a 
radial spherically symmetric average will give an average Hubble constant whose value decreases from the maximum 
at the 30/z -1 Mpc scale to the dressed global average value, as the scale of homogeneity is approached at roughly the 
baryon acoustic oscillation (BAO) scale of HO/i Mpc. This predicted effect could account for the Hubble bubble 
UJ1 and more detailed studies of the scale dependence of the local Hubble flow lfl2tl . 

In fact, the variance of the local Hubble flow below the scale of homogeneity should correlate strongly to observed 
structures in a manner which has no equivalent prediction in FLRW models. 



FUTURE OBSERVATIONAL TESTS 

There are two types of potential cosmological tests that can be developed; those relating to scales below that of 
statistical homogeneity as discussed above, and those that relate to averages on our past light cone on scales much 
greater than the scale of statistical homogeneity. The second class of tests includes equivalents to all the standard 
cosmological tests of the standard FLRW model with Newtonian perturbations. This second class of tests can be further 
divided into tests which just deal with the bulk cosmological averages (luminosity and angular diameter distances etc), 
and those that deal with the variance from the growth of structures (late epoch integrated Sachs- Wolfe effect, cosmic 
shear, redshift space distortions etc). Here I will concentrate solely on the simplest tests which are directly related to 
luminosity and angular diameter distance measures. 

In the timescape cosmology we have an effective dressed luminosity distance 

d L =a (l+z)r w , (6) 

where a Q = f Q l a Q , and 

r w = r(i-/v) 1/3 /°^ )(1 _ * (f , ))1/3fl - (f , r (7) 

We can also define an effective angular diameter distance, d A , and an effective comoving distance, D, to a redshift z in 
the standard fashion 

D d, 

A = T+z = (T+5]2 ' (8) 
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FIGURE 1. The effective comoving distance H Q D(z) is plotted for the best-fit timescape (TS) model, with f V Q = 0.762, (solid 
line); and for various spatially flat ACDM models (dashed lines). The parameters for the dashed lines are (i) £2 M() = 0.249 (best-fit 
to WMAP5 only QJ]); (ii) fl M) = 0.279 (joint best-fit to Snela, BAO and WMAP5); (iii) £l m = 0.34 (best-fit to Riess07 Snela 
only [14]). Panel (a) shows the redshift range z < 6, with an inset for z < 1.5, which is the range tested by current Snela data. Panel 
(b) shows the range z < 1100 up to the surface of last scattering, tested by WMAP5. 

A direct method of comparing the distance measures with those of homogeneous models with dark energy, is to 
observe that for a standard spatially flat cosmology with dark energy obeying an equation of state P D = w(z)p D , the 
quantity 



does not depend on the value of the Hubble constant, H Q , but only directly on £2 MQ = 1 — ii D0 - Since the best-fit values 
of Hq are potentially different for the different scenarios, a comparison of H Q D curves as a function of redshift for the 
timescape model versus the ACDM model gives a good indication of where the largest differences can be expected, 
independently of the value of H Q . Such a comparison is made in Fig.Q] 

We see that as redshift increases the timescape model interpolates between ACDM models with different values of 
fi MQ . For redshifts z < 1.5 D TS is very close to D ACBM for the parameter values (£2 MQ ,fi A0 ) = (0.34,0.66) (model 
(iii)) which best-fit the Riess07 supernovae (Snela) data [ 14] only, by our own analysis. For very large redshifts 
that approach the surface of last scattering, z < 1100, on the other hand, D TS very closely matches D ACBM for 
the parameter values (Q. m ,D. AQ ) = (0.249,0.751) (model (i)) which best-fit WMAP5 only 0]. Over redshifts 
2 < z < 10, at which scales independent tests are conceivable, Z) Tg makes a transition over corresponding curves 
of D ACDM with intermediate values of (£1 M0 ,£1 A0 ). The D ACDM curve for joint best-fit parameters to Snela, BAO 
measurements and WMAP5 11311 . (n M0 ,£2 A0 ) = (0.279,0.721) is best-matched over the range 5 < z < 6, for example. 

The difference of D TS from any single D ACDU curve is perhaps most pronounced in the range 2 < z < 6, which may 
be an optimal regime to probe in future experiments. Gamma-ray bursters (GRBs) now probe distances to redshifts 
z < 8.3, and could be very useful. A considerable amount work of work has already been done on Hubble diagrams for 
GRBs. (See, e.g., B15I1 .) Much more work is needed to nail down systematic uncertainties, but GRBs may eventually 
provide a definitive test in future. An analysis of the timescape model Hubble diagram using 69 GRBs has just been 
performed by Schaefer flrjll . who finds that it fits the data better than the concordance ACDM model, but not yet by a 
huge margin. As more data is accumulated, it should become possible to distinguish the models if the issues with the 
standardization of GRBs can be ironed out. 



It should be noted that the shape of the H Q D curves depicted in Fig. Q] represent the observable quantity one is 
actually measuring when some researchers loosely talk about "measuring the equation of state". For spatially flat 
dark energy models, with H Q D given by (O, one finds that the function w(z) appearing in the fluid equation of state 




(9) 



The effective "equation of state' 




FIGURE 2. The artificial equivalent of an equation of state constructed using the effective comoving distance d 1 Ot . plotted for 
the timescape tracker solution with best-fit value / v n = 0.762, and two different values of S1 M q- (a) the canonical dressed value 

a ito = 2(1 -/vo)(2+/vo) = 0.33; (b) £l m = 0.279. 



P D = w{z)p D is related to the first and second derivatives of (O by 

j(l+;)P'-'P" + l 

"M-iWi+^-i <10) 

where prime denotes a derivative with respect to z. Such a relation can be applied to observed distance measurements, 
regardless of whether the underlying cosmology has dark energy or not. Since it involves first and second derivatives 
of the observed quantities, it is actually much more difficult to determine observationally than directly fitting H D(z). 

The equivalent of the "equation of state", w(z), for the timescape model is plotted in Fig. [2] The fact that w(z) 
is undefined at a particular redshift and changes sign through ±°° simply reflects the fact that in ( TTOb we are 
dividing by a quantity which goes to zero for the timescape model, even though the underlying curve of Fig. Q~|is 
smooth. Since one is not dealing with a dark energy fluid in the present case, w (z) simply has no physical meaning. 
Nonetheless, phenomenologically the results do agree with the usual inferences about w (z) for fits of standard dark 
energy cosmologies to Snela data. For the canonical model of Fig. [2ja) one finds that the average value of w(z) — — 1 
on the range z < 0.7, while the average value of w(z) < — 1 if the range of redshifts is extended to higher values. The 
w = — 1 "phantom divide" is crossed at z — 0.46 for / v q ~ 0.76. One recent study [ 17] finds mild 95% evidence for 
an equation of state that crosses the phantom divide from w > — 1 to w < - lin the range 0.25 < z < 0.75 in accord 
with the timescape expectation. By contrast, another study 111 80 at redshifts z < 1 draws different conclusions about 
dynamical dark energy, but for the given uncertainties in w(z) the data is consistent with Fig. [TJa) as well as with a 
cosmological constant |Q]]. 

The fact that w(z) is a different sign to the dark energy case for z > 2 is another way of viewing our statement above 
that the redshift range 2 < z < 6 may be optimal for discriminating model differences. 



The H(z) measure 

Further observational diagnostics can be devised if the expansion rate H(z) can be observationally determined as 
a function of redshift. Recently such a determination of H(z) at z = 0.24 and z = 0.43 has been made using redshift 
space distortions of the BAO scale in the ACDM model IU9I1 . This technique is of course model dependent, and the 
Kaiser effect would have to be re-examined in the timescape model before a direct comparison of observational results 
could be made. A model-independent measure of H(z), the redshift time drift test, is discussed below. 

In Fig.|3]we compare H(z)/H for the timescape model to spatially flat ACDM models with the same parameters 
chosen in Fig. Q] The most notable feature is that the slope of H(z)/H Q is less than in the ACDM case, as is to be 
expected for a model whose (dressed) deceleration parameter varies more slowly than for ACDM. 
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FIGURE 3. The function H~ x H(z) for the timescape model with f V Q = 0.762 (solid line) is compared to H~ l H(z) for three 
spatially flat ACDM models with the same values of (fi MQ , fi AQ ) as in Fig.[T](dashed lines). 



The Om(z) measure 

Recently a number of authors Ir2(il,l2"ill22tl have discussed various roughly equivalent diagnostics of dark energy. For 
example, Sahni, Shafieloo and Starobinsky |2j] , have proposed a diagnostic function 



Om{£) = 



H 2 (z) 



- 1 



(11) 



on account of the fact that it is equal to the constant present epoch matter density parameter, £i M0 , at all red- 
shifts for a spatially flat FLRW model with pressureless dust and a cosmological constant. However, it is not con- 
stant if the cosmological constant is replaced by other forms of dark energy. For general FLRW models, H(z) = 
i/i r> H-nJi-i .... - ..;„„!.■, ,w;,..,.;,..-,. «r ni-l Thus the diagnostic dlU is easier to 



[Z)'(z)] 'y'l + Q, kQ H 2 D 2 (z), which only involves a single derivatives of D(z) 
reconstruct observationally than the equation of state parameter, w(z). 




(a) 

FIGURE 4. The dark energy diagnostic Om(z) of Sahni, Shafieloo and Starobinsky 12111 plotted for the timescape tracker solution 
with best-fit value / v o = 0.762 (solid line), and 1(7 limits (dashed lines) from ref. Ildl : (a) for the redshift range < z < 1.6 as 
shown in ref. [23]; (b) for the redshift range < z < 6. 

The quantity Om(z) is readily calculated for the timescape model, and the result is displayed in Fig. [4] What is 
striking about Fig. [4] as compared to the curves for quintessence and phantom dark energy models as plotted in ref. 
12111. is that the initial value 



Om(0) 



2 n'\ 

3 H lo 



2(8/ v 3 -3/ v 2 +4)(2+/ v0 ) 



(4/ v 2 +/v0+4) 2 



(12) 



is substantially larger than in the spatially flat dark energy models. Furthermore, for the timescape model Om(z) does 
not asymptote to the dressed density parameter i2 M0 in any redshift range. For quintessence models Om(z) > & M0 , 
while for phantom models Om(z) < Q MQ , and in both cases Om(z) — > £2 M0 as z — > °°. In the timescape model, 
Om(z) > £2 M) — 0.33 for z < 1 .7, while Om(z) < i2 M0 for z > 1 .7. It thus behaves more like a quintessence model for 
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FIGURE 5. (a) The Alcock-Paczynski test function / Ap = HD/r, and (b) the BAO radial test function H Q D V = H Q Df A ^. In 
each case the timescape model with / v o = 0.762 (solid line) is compared to three spatially flat ACDM models with the same values 
of (£l M0 ,Q, AQ ) as in Fig. [T] (dashed lines). 



low z, in accordance with Fig. [2] However, the steeper slope and the different large z behaviour mean the diagnostic is 
generally very different to that of typical dark energy models. For large z, & MQ < Om(°°) < £2 M0 , if / v o > 0.25. 

Interestingly enough, a recent analysis of Snela, BAO and CMB data l23ll for dark energy models with two different 
empirical fitting functions for w(z) gives an intercept Om(0) which is larger than expected for typical quintessence or 
phantom energy models, and in the better fit of the two models the intercept (see Fig. 3 of ref. [23]) is close to the value 
expected for the timescape model, which is tightly constrained to the range 0.638 < Om(0) < 0.646 if / v o = 0.76^Q Qg. 



The Alcock-Paczynski test and baryon acoustic oscillations 



Some time ago Alcock and Paczyhski devised a test H24I1 which relies on comparing the radial and transverse proper 
length scales of spherical standard volumes comoving with the Hubble flow. This test, which determines the function 



AP 



89_ 

8z 



HD 



(13) 



was originally conceived to distinguish FLRW models with a cosmological constant from those without a A term. 
The test is free from many evolutionary effects, but relies on one being able to remove systematic distortions due to 
peculiar velocities. 

Current detections of the BAO scale in galaxy clustering statistics IT251 l26ll can in fact be viewed as a variant of 
the Alcock-Paczynski test, as they make use of both the transverse and radial dilations of the fiducial comoving BAO 
scale to present a measure 



zD z 
H(z) 



1/3 



Df, 



-1/3 
AP " 



(14) 



In Fig.[5]the Alcock-Paczynski test function ( fT3l and BAO scale measure ( TBI of the timescape model are compared 
to those of the spatially flat ACDM model with different values of (i2 A0 ,£l A0 ). Over the range of redshifts z < 1 studied 
currently with galaxy clustering statistics, the / Ap curve distinguishes the timescape model from the ACDM models 
much more strongly than the D y test function. In particular, the timescape / Ap has a distinctly different shape to that 
of the ACDM model, being convex. The primary reason for use of the integral measure (TBI has been a lack of data. 
Future measurements with enough data to separate the radial and angular BAO scales are a potentially powerful way 
of distinguishing the timescape model from ACDM. 

Recently Gaztanaga, Cabre and Hui [ 19] have made the first efforts to separate the radial and angular BAO scales 
in different redshift slices. Although they have not yet published separate values for the radial and angular scales, 
their results are interesting when compared to the expectations of the timescape model. Their study yields best-fit 
values of the present total matter and baryonic matter density parameters, i2 M0 and £2g , which are in tension with 
WMAP5 parameters fit to the ACDM model. In particular, the ratio of non-baryonic cold dark matter to baryonic 
matter has a best-fit value i2 r „/i2„ = (Q. M() Q. R() of 3.7 in the 0.15 < z < 0.3 sample, 2.6 in the 0.4 < z <0.47 



FIGURE 6. Left panel: The (in)homogeneity test 
the timescape tracker solution with best-fit value 
the equivalent curves SS = Q, k() (H () D) 2 for two 
(a) a m = 0.28, fi A0 = 0.71, Cl k0 = 0.01; (b) Q M0 = 0.28, £2 A0 = 0. 
Right panel: The (in)homogeneity test function ff(z) is plotted for tl 



function 3§(z) = [HD 1 ] 2 - 1 is plotted for 
fvO = 0.762 (solid line), and compared to 
different ACDM models with small curvature: 
73, a k0 = -0.01. 
le fvO — 0.762 tracker solution. 



sample, and 3.6 in the whole sample, as compared to the expected value of 6.1 from WMAP5. The analysis of the 
3-point correlation function yields similar conclusions, with a best fit 12711 fi M0 = 0.28 ±0.05, Q. B0 = 0.079 ±0.025. 
By comparison, the parameter fit to the timescape model of ref. IToll yields dressed parameters Q M0 = 0.33^g'[g, 

£l B0 = O.O8OI0 013, an d a ra ti° ^co/^bo = 31^2 4- Since other forms of dark energy are not generally expected to 
give rise to a renormalization of the ratio of non-baryonic to baryonic matter, this is encouraging for the timescape 
model. 



Test of (in)homogeneity 

Recently Clarkson, Bassett and Lu If28ll have constructed what they call a "test of the Copernican principle" based 
on the observation that for homogeneous, isotropic models which obey the Friedmann equation, the present epoch 
curvature parameter, a constant, may be written as 

a a JJ&Mfl (15) 

for all z, irrespective of the dark energy model or any other model parameters. Consequently, taking a further derivative, 
the quantity 

<*?(z) = 1 +H 2 (DD" -D' 2 ) + HH'DD 1 (16) 

must be zero for all redshifts for any FLRW geometry. 

A deviation of ^(z) from zero, or of (T5[ from a constant value, would therefore mean that the assumption of 
homogeneity is violated. Although this only constitutes a test of the assumption of the Friedmann equation, i.e., of the 
Cosmological Principle rather than the broader Copernican Principle adopted in ref. 01, the average inhomogeneity 
will give a clear and distinct prediction of a non-zero ^(z) for the timescape model. 

The functions (15[ and ( fT6b are computed in ref. [1]. Observationally it is more feasible to fit (T5[ which involves one 
derivative less of redshift. In Fig.|6]we exhibit both ^(z), and also the function 3§(z) = [HD 1 ] 2 — 1 from the numerator 
of ( TT3T > for the timescape model, as compared to two ACDM models with a small amount of spatial curvature. A 
spatially flat FLRW model would have 5§{z) = 0. In other FLRW cases £%(z) is always a monotonic function whose 
sign is determined by that of £2 W) - An open A = universe with the same £2 M) would have a monotonic function 38{z) 
very much greater than that of the timescape model. 



Time drift of cosmological redshifts 



For the purpose of the Om(z) and (in)homogeneity tests considered in the last section, H (z) must be observationally 
determined, and this is difficult to achieve in a model independent way. There is one way of achieving this, however, 



namely by measuring the time variation of the redshifts of different sources over a sufficiently long time interval l29tl . 
as has been discussed recently by Uzan, Clarkson and Ellis [30]. Although the measurement is extremely challenging, 
it may be feasible over a 20 year period by precision measurements of the Lyman-a forest in the redshift range 
2 < z < 5 with the next generation of Extremely Large Telescopes i3lll . 

In ref . IH an analytic expression for H~ 1 j| is determined, the derivative being with respect to wall time for observers 
in galaxies. The resulting function is displayed in Fig. |7]for the best-fit timescape model with / v o = 0.762, where it is 
compared to the equivalent function for three different spatially flat ACDM models. What is notable is that the curve 
for the timescape model is considerably flatter than those of the ACDM models. This may be understood to arise from 
the fact that the magnitude of the apparent acceleration is considerably smaller in the timescape model, as compared to 
the magnitude of the acceleration in ACDM models. For models in which there is no apparent acceleration whatsoever, 
one finds that H Q 1 ^ is always negative. If there is cosmic acceleration, real or apparent, at late epochs then H~ l ^ 
will become positive at low redshifts, though at a somewhat larger redshift than that at which acceleration is deemed 
to have begun. 
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FIGURE 7. The function H Q 1 j| for the timescape model with / v q = 0.762 (solid line) is compared to 1 j| for three spatially 
flat ACDM models with the same values of (£2 MQ) i2. ) as in Fig.[TJ(dashed lines). 

Fig.[7]demonstrates that a very clear signal of differences in the redshift time drift between the timescape model and 
ACDM models might be determined at low redshifts when H~ l || should be positive. In particular, the magnitude of 

3§ i s considerably smaller for the timescape model as compared to ACDM models. Observationally, however, it 
is expected that measurements will be best determined for sources in the Lyman a forest in the range, 2 < z < 5. At 
such redshifts the magnitude of the drift is somewhat more pronounced in the case of the ACDM models. For a source 
at z = 4, over a period of 8% = 10 years we would have Sz — — 3.3 x 10~ 10 for the timescape model with / v o = 0.762 
and Hq = 61.7 km sec -1 Mpc~'. By comparison, for a spatially flat ACDM model with H Q = 70.5 km sec -1 Mpc~' a 
source at z = 4 would over ten years give 8z= -4.7 x 10~ 10 for (Q. m ,Q. AQ ) = (0.249,0.751), and <5z = -7.0 x 10~ 10 
for (n M0 ,n A0 ) = (0.279,0.721). 



DISCUSSION 

The tests outlined here demonstrate several lines of investigation to distinguish the timescape model from models of 
homogeneous dark energy. The (in)homogeneity test of Bassett, Clarkson and Lu is definitive, since it tests the validity 
of the Friedmann equation directly. 

In performing these tests, however, one must be very careful to ensure that data has not been reduced with built-in 
assumptions that use the Friedmann equation. For example, current estimates of the BAO scale such as that of Percival 
et al. [26] do not determine D y directly from redshift and angular diameter measures, but first perform a Fourier space 
transformation to a power spectrum, assuming a FLRW cosmology. Redoing such analyses for the timescape model 
may involve a recalibration of relevant transfer functions. 

In the case of supernovae, one must also take care as compilations such as the Union [32] and Constitution ll33ll 
datasets use the SALT method to calibrate light curves. In this approach empirical light curve parameters and cosmo- 
logical parameters - assuming the Friedmann equation - are simultaneously fit by analytic marginalisation before the 
raw apparent magnitudes are recalibrated. As Hicken et al. discuss [HH, a number of systematic discrepancies exist 
between data reduced by the SALT, SALT2, MLCS31 and MLCS17 techniques even within the ACDM model. In 



the case of the timescape model, we find considerable differences between the different approaches ll34ll . In principle, 
at present there appear to be enough supernovae to decide between the ACDM and timescape models on Bayesian 
evidence, but one is led to different conclusions depending on how the data is reduced. It is therefore important that 
the systematic issues are unravelled. 

The value of the dressed Hubble constant is also an observable quantity of considerable interest. A recent determi- 
nation of H by Riess et al. IB5I1 poses a challenge for the timescape model. However, it is a feature of the timescape 
model that a 17-22% variance in the apparent Hubble flow will exist on local scales below the scale of statistical ho- 
mogeneity, and this may potentially complicate calibration of the cosmic distance ladder. Further quantification of the 
variance in the apparent Hubble flow in relationship to local cosmic structures would provide an interesting possibility 
for tests of the timescape cosmology for which there are no counterparts in the standard cosmology. 
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